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Abstract. The objective of this effort is to perform design analyses for a non-nuclear hot-hydrogen materials tester, as a first 
step towards developing efficient and accurate multiphysics, thermo-fluid computational methodology to predict 
environments for hypothetical solid-core, nuclear thermal engine thrust chamber design and analysis. The computational 
methodology is based on a multidimensional, finite-volume, turbulent, chemically reacting, thermally radiating, unstructured- 
grid, and pressure-based formulation. The multiphysics invoked in this study include hydrogen dissociation kinetics and 
thermodynamics, turbulent flow, convective, and thermal radiative heat transfers. The goals of the design analyses are to 
maintain maximum hot-hydrogen jet impingement energy and to minimize chamber wall heating. The results of analyses on 
three test fixture configurations and the rationale for final selection are presented. The interrogation of physics revealed that 
reactions of hydrogen dissociation and recombination are highly correlated with local temperature and are necessary for 
accurate prediction of the hot-hydrogen jet temperature. 
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INTRODUCTION 

Nuclear thermal propulsion (NTP) may open up the solar system to far broader and faster exploration than is now 
possible with chemical propulsion. The feasibility of NTP systems was established by extensive testing in the 
Rover/NERVA programs and the technical merits of NTP have been identified in numerous studies (Koenig, 1986). In 
a NTP system, the propellant is heated by a nuclear reactor instead of combustion. The best NTP propellant is hydrogen 
with which high exhaust speed can be achieved with its low molecular weight, resulting in high thrust per unit mass of 
propellant consumed. Under the operating temperature inside a nuclear reactor, often at 2,500-3,000 K, the heated 
hydrogen further dissociates into atomic hydrogen - at a formula weight only half that of molecular hydrogen, indicating 
an even greater thrust may be realized. The heat transfer efficiency and degree of hydrogen dissociation are therefore 
two important performance factors for NTP. On the other hand, the need to push fuel element temperature to extremes in 
order to maximize performance also intensifies hydrogen induced corrosion rates, which are known to increase in direct 
proportion to reactor operating temperature (Lyon, 1973). In order to develop candidate high temperature fuel materials 
that would be compatible with the hot-hydrogen environment of a high performance solid core NTP engine, a non- 
nuclear test effort entitled “Hot Hydrogen Materials and Component Development” is underway at NASA’s Marshall 
Space Flight Center (MSFC), while a parallel task entitled “Multiphysics Thrust Chamber Modeling” is also in progress 
in order to develop a computational methodology capable of predicting the thermal-fluid environment in a nuclear 
thermal engine thrust chamber. This computational methodology is based on an Unstructured-grid Navier-Stokes 
Intemai-extemai computational fluid dynamics (CFD) Code (UNIC). Physical and numerical models pertinent to solid- 
core NTP will be developed and implemented. The hot gas and material temperatures to be measured in the hot- 
hydrogen materials development tests will be used to benchmark the UNIC code. In this effort, the design part of the 
UNIC capability was demonstrated by performing design analyses for three hot hydrogen materials tester configurations 
and a final configuration was recommended. The effect of configuration changes on design goals was presented and the 
effects of hydrogen dissociation and thermal radiation on the computed thermophysics were investigated. 



TEST FIXTURE BASELINE DESIGN 


The baseline design of the apparatus for testing tubular 
fuel materials is shown in Fig. 1. It is directly mated to 
an arc-heater (not shown) that provides hot-hydrogen 
flow. Optical ports are fitted to allow real-time 
pyrometer and laser diagnostics measurements for 
material surface temperature and centerline gas 
temperature, respectively. As the hot-hydrogen jet 
travels the length of the chamber, it loses energy to the 
water-cooled copper chamber and sample material. As 
the sample material heats up, it loses energy to the 
colder chamber wall by thermal radiation. The first 
design goal is therefore minimizing the energy loses of 
the hydrogen jet. As the temperature of the hydrogen jet 
reaches that of NTP operating condition, the heat load 
imposed on the copper chamber may exceed its design 
limit, due to both the convective heating from hydrogen 
jet and bombardment of thermal radiation from hot 
sample. The second design goal is therefore to minimize 
the heat load on the copper chamber wall. 



FIGURE 1. Baseline Design of the Test Fixture. 


COMPUTATIONAL METHODOLOGY 


The design analyses were performed using UNIC CFD methodology. UNIC is a multidimensional, finite-volume, 
turbulent, chemically reacting, thermally radiating, unstructured-grid, and pressure-based formulation. The all-speed flow nature 
makes it suitable for describing the thermophysics encountered in a NTP environment. UNIC was developed recently through 
several activities, namely the launch vehicle base-heating (Chen, 2001), Laser propulsion (Wang, 2002), and stage- 
separation (Chen, 2002). Most recently, UNIC was benchmarked for steady-state heat transfer and thrust 
performance for the Space Shuttle Main Engine (SSME) (Wang, 2004) and SSME transient thrust chamber side load 
(Wang, 2005). To describe the thermal-fluid design environment of the tester, UNIC solves the time-varying 
transport equations of continuity, species continuity, momentum, total enthalpy, turbulent kinetic energy, and 
turbulent kinetic energy dissipation: 
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A predictor and corrector solution algorithm was employed to provide coupling of the fluid governing equations. A 
second-order central-difference scheme was employed to discretize the diffusion fluxes and source terms of the 
governing equations. For the convective terms, a second-order upwind total variation diminishing difference scheme 
was used in this effort. A point-implicit (operator splitting) method was employed to solve the finite-rate chemistry 
system. An extended k-e turbulence model (Chen, 1987) was used to describe the turbulence. 

The convective heat transfer follows the modified Newtonian law 

Q c =(pu r IT + )(h w -h p - R(u 2 p l2 )) (7) 

The radiative heat transfer is analyzed by solving the radiative transfer equation 

(Q * V)/(r, Q) = - Kl(r , Q) + Kl b (r) (8) 

Discrete ordinate method was used to solve the radiative transfer equation. The radiative heat flux is given by the 
integration of the wall leaving radiative intensities 

Q r = J/(r,Q’)|«U’|dQ“ (9) 

n .^“<0 

Details of the numerical algorithm can be found in (Wang, 2003 and 2004). 

Boundary Conditions 

No-slip condition was applied to the solid walls. Fixed mass flow rate boundary condition was used at the inlet, and 
mass conservation boundary condition was used at the exit. For a conservative calculation of the chamber wall heat flux, 
a fixed temperature of 400 K was estimated for the chamber wall and the hot-hydrogen temperature was set at 3500 K. 
For the shield configuration, the wall temperature of the upper shield facing the chamber wall was estimated to be 2600 
K from a separate one-dimensional heat transfer calculation. Adiabatic condition was applied to walls of sample, funnel 
(for funnel configuration), lower shield and flanges, and graphite 
insert. The hydrogen inlet mass flow rate was 10 g/s and the chamber 
pressure was 35 atm. An emissivity of 0.6 was applied to the copper 
chamber wall, whereas an emissivity of 0.4 was applied to the shield. 

The emissivity for the rest of the solid walls was set at 0.9. The inlet 
boundary was considered as a radiating wall to approximate the 
radiation from the arc-heater section. A series of pre-calculations were 
performed on the baseline case to iterate the inlet temperature and 
species concentrations such that the inlet species concentrations 
correspond to a state of temperature at 3500 K. 

DESIGN ANALYSIS 

Multiphysics invoked in this effort include turbulence, hydrogen 
reactions (Wang, 2001) and its associated thermodynamics, 
convective heat transfer, and surface radiation. The goals of this 
design are to maintain the 3500 K hot-hydrogen jet as long as possible 
as it impinges on the sample and to hold the overall chamber wall heat 
flux below a pre-determined 5.0x1 0 6 W/m 2 . 

Baseline Design 




The hot hydrogen (H 2 ) environment of the baseline design was computed and Figure 2 shows the computed temperature, 
hydrogen atom (H) mass fraction contours and flow streamlines. Temperature contours show the temperature drop of the 
hot-hydrogen jet, as it loses energy to the colder copper chamber wall due to both convection and sample surface 
radiation processes. As the hot jet cools, H recombines to become H 2 and it can be seen from the H contours that higher 
concentration of H is correlated with hotter local temperature. The streamline plot shows an expanding hot-hydrogen jet 
impinging on and flowing around the sample block, and later exhausting into the convergent exit section. The sample 
block is heated by the hot-hydrogen jet and its integrity under such a high temperature environment is the subject of the 
test. A large recirculation zone appears in the divergent section of the chamber, while a small recirculation region forms 
behind the sample holder; both of which are strongly affected by the turbulence. Note although the plots in Fig. 1 
showing a distinctive strong hot-hydrogen jet and two recirculated flow regions, the pressure and Mach number contours 
are fairly uniform inside the chamber (not shown) for this design, due to the largely low subsonic flow field. 


Figure 3 shows the computed centerline temperature profiles from inlet to the sample front face for the baseline 
design. The end of the curves is the stagnation point. The results of four cases are presented in here, with combinations 
of frozen and finite-rate chemistries and with and without surface radiation. Finite-rate chemistry considers the H 2 
dissociation and H recombination, while the frozen chemistry freezes the inlet flow composition. The centerline line 
temperatures appear to maintain its value until the energy dissipation starts to take over and they drop rapidly until 
reaching the stagnation point. Without surface radiation, the central surface (stagnation) temperatures for frozen and 
finite-rate chemistries are 2954.1 and 3170.8 K, respectively. The temperature computed by frozen chemistry is 216.7 K 
lower because thermodynamic properties were computed based on an unphysical frozen composition, indicating the 
importance of finite-rate chemistry in arriving accurate thermal-fluid solution. With surface radiation, the central sample 
surface temperatures drop further due to surface energy exchange with the colder chamber wall, resulting in stagnation 
temperatures of 2805.6 and 3052.8 K for frozen and finite-rate chemistries, respectively. Figure 3 demonstrates the 
importance of the inclusion of both hydrogen dissociation and surface radiation in the design analysis. Figure 3 also 
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FIGURE 3. Computed Centerline Temperature FIGURE 4. Computed Total, Convective, and 

Profiles from Inlet to Stagnation Point for Baseline Radiative Heat Fluxes on Chamber Wall for Baseline 

Design Cases. Design. 

shows that the approximately 450 K temperature drop from the inlet to the stagnation point indicating too much energy 
dissipation and the baseline design needs to be improved. In addition, the 118.0 K temperature drop with and without 
radiation indicates a radiation shield may be in order. 


total, convective, and radiative heat fluxes along the chamber wall between anode and 


graphite insert sections. Although the surface radiation caused 1 18.0 K drop in stagnation temperature, the radiative heat 
flux is much lower than that of the convective heat flux in comparison. This is because the area of the sample radiating 
surface is much smaller than that of the fluid-chamber contact surface. And neither H 2 nor H is participating medium for 
radiation. Figure 4 also reveals that the total heat flux is mostly above 5xl0 6 W/m 2 along the chamber surface, with a 
peak heat flux at about lxlO 7 W/m 2 at the hot-hydrogen flow impingement point on the chamber, indicating the baseline 




design may have overheating issue as well. 


Funnel Design 


The second design involves placing a funnel in front of the coupon. A funnel would concentrate the hot-hydrogen flow 
and possibly shield part of the convective energy exchange between the hot jet and the cooled chamber, while also serve 
as a partial radiation-shield. 


Figure 5 shows the computed temperature, H concentration and 
streamlines for the funnel design. It can be seen that the 
expansion of the hot-hydrogen jet in the divergent section is 
forced to converge into a more energy concentrated jet, by the 
funnel. This is evident by the slower decreasing slope of the 
centerline temperature profiles in Fig. 6, comparing to the steep 
decreasing slope of those in Fig. 3. The stagnation temperatures 
are 3172.0 K (with radiation) and 3191.7 K (without radiation). 
The stagnation temperature drop with and without radiation is 
now only 19.7 K, indicating that the funnel indeed serves as a 
partial radiation shield. And again the higher concentration of 
H associates with higher temperatures. 

The flowfield of the funnel design however, is not as uniform 
as that of the baseline design. As the higher energy, 
concentrated hot-hydrogen jet impinges on the sample, instead 
of flowing around it smoothly as in the baseline case, it forms a 
stronger jet impinging on the chamber wall, as indicated in the 
streamline plot in Fig. 5. That results in a peak total chamber 
wall heat flux of 1.54xl0 7 W/m 2 , exceeding the maximum heat 
flux desirable, as indicated in Fig. 7. In fact, the overall total 
heat flux exceeds the desired 5xl0 6 W/m 2 , except at the root of 
the funnel. 



FIGURE 5. Computed Temperature, H Mass Fraction 
Contours and Streamlines for Funnel Design. 




FIGURE 6. Computed Centerline Temperature Profiles FIGURE 7. Computed Total, Convective, and Radiative 

from Inlet to Sample Front Face for Funnel Design Cases. Heat F^xes on Chamber Wall for Funnel Design. 



Shield Design 


The third design involves a Tungsten shield that almost completely shields the chamber wall from the convective heating 
and is supported by a 360° Flange. The beginning of the shield is closed off at the intersection of the chamber and the 


anode cap so that no hot flow can get in between the shield 
and chamber wall Without the 360° Flange and the close- 
off, a separate calculation indicated that a flow pumping 
effect is invoked and the hot hydrogen is sucked into the 
region between the chamber the shield. 

Figure 8 shows the computed temperature, H 
concentration contours and streamlines for the shield 
design. It can be seen from the temperature contours that 
the flow is now uniformly hotter than that of the baseline 
design, because the convective heat loss is mostly 
shielded. The H jet is now propagating the longest while 
wrapping around the coupon, indicating minimum energy 
dissipation for the hot jet. The H concentration between 
the chamber and the shield is low because of the low local 
temperature. The flow streamlines look very similar to 
those of the baseline design except the flow speed is a little 
faster due to the smaller effective flowing area, resulting in 
a smaller recirculation bubble at the divergent section. 

Figure 9 shows that the dropping rate of the centerline 
temperatures is now very slow, comparing to those of Fig. 
3 and 6. The stagnation temperatures are 3426.0 K with 
radiation, and 3440.7 K without radiation. The computed 
stagnation temperatures with and without thermal radiation 
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FIGURE 8. Computed Temperature, H Mass Fraction 
Contours and Streamlines for Shield Design. 
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FIGURE 9 Computed Centerline Temperature Profiles 
for Shield Configuration Cases. 


FIGURE 10. Computed Total, Convective, and Radiative 
Heat Fluxes on Chamber W all for Shield Design. 


calculation for all three designs are summarized in Table 1. The smallest stagnation temperature drop of 14.7 K with and 
without radiation indicating this shield is also a very effective radiation shield. Figure 10 shows the average total heat 
flux is now an acceptable 2.5xl0 6 W/m 2 , lower than the 5xl0 6 W/m 2 limit. There is a peak total heat flux of 5.65xl0 6 




W/m 2 right behind the close-off shield that is probably caused by the pre-determined 2600 K for the chamber-facing 
shield surface. Since the close-off part of the shield is connected to the cold chamber wall at 400 K, it is believed that the 
shield surface temperature at the close-off should be 
lower due to conduction, resulting in lower local gas 
temperature, and therefore lower peak heat flux if a 
coupled conjugate heat transfer calculation was 
conducted. The streamlines plot show a small 
recirculation bubble at the end of the shield, resulting in 
a second peak flux of 3.46xl0 6 W/m 2 , also lower than 
the threshold. The shield configuration is therefore the 
design satisfying both design criteria. Figure 1 1 shows 
the test fixture for testing rod-shaped fuel materials with 
this shield design. 


Table 1. Computed Stagnation Temperatures for Three Test 


Design 

Without 
Radiation, K 

With 

Radiation, K 

A, K 

baseline 

3170.8 

3052.8 

118.0 

funnel 

3191.7 

3172.0 

19.7 

shield 

3440.7 

3426.0 

14.7 



FIGURE 11. Shield Design. 


CONCLUSION 

Multiphysics computational analyses were performed to support the design of test fixtures for non-nuclear testing of fuel 
materials. Multiphysics invoked include hydrogen dissociation kinetics and thermodynamics, turbulent flow, convective, 
and thermal radiative heat transfers. The design criteria were to maximize the hot-hydrogen temperature while impinging 
on the test coupon, and to minimize the total heat flux imposed on the copper chamber wall. Results of the design 
analyses indicate that the shield design provides the lowest chamber wall total heat flux and the highest hot hydrogen 
delivering temperature, thereby is superior to the baseline and funnel designs. The effects of sample surface radiation on 
chamber heat flux and sample stagnation temperature were presented. It was also found that hydrogen dissociation and 
recombination reactions are highly correlated with local temperature and are very important in accurate prediction of the 
chamber thermo-flowfield. 


NOMENCLATURE 

C h C 2 ,C 3 ,C^ turbulence modeling constants, 1.15, 1.9, 0.25, and 0.09. 
C p - heat capacity (m 2 /(s 2 K)) 

D - diffusivity (m 2 /s) 

H = total enthalpy (J) 

h = static enthalpy (J) 

/ = radiative intensity (W/(m 2 sr)) 

K - thermal conductivity (W/(m K)) 

k = turbulent kinetic energy (m 2 /s 2 ) 

p = pressure (atm) 

Q = heat flux (W/m 2 ) 

R - recovery factor 

r = location coordinate (m) 

T = temperature (K) 

F* = non-dimensional temperature 

t = time (s) 

ui - mean velocities in three directions (m/s) 

u x - non-dimensional wall friction velocity 


* » 


x = Cartesian coordinates (m) 

a = species mass fraction 

8 = turbulent kinetic energy dissipation rate (m 2 /s 3 ) 

6 = energy dissipation contribution (kg/(m 2 s 3 ) 

k = absorption coefficient (atm m)' 1 

ft = viscosity (kg/(m s)) 

p t - turbulent eddy viscosity (=pC^ 2 /e) (kg/(m s)) 

n - turbulent kinetic energy production (m 2 /s 3 ) 

p = density (kg/m 3 ) 

a = turbulence modeling constants 

r = shear stress (kg/(m s 2 ) 

Q = direction vector. Cl denotes the leaving radiative intensity direction (sr) 

co = chemical species production rate (kg/(m 3 s)) 


Subscripts 


b 

c 

cl 

r 

P 

t 

w 


= black body 
= convective 
= centerline 
= radiative 
= off-wall point 
= turbulent flow 
= wall point 
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Governing Equations - continued 
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NASA MSFC Multi-Gas Arc-Heater Facility 
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Hot Hydrogen Materials and Component Tester 
Design Analysis - Baseline Design 
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Design Analysis - Baseline Design 
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Hot Hydrogen Materials and Component Tester 
Design Analysis - Baseline Design 
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Hot Hydrogen Materials and Component Tester 
Design Analysis - Funnel Design 
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Hot Hydrogen Materials and Component Tester 
Design Analysis - Shield Design 
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Hot Hydrogen Materials and Component Tester 

Desian Analvsis 



Shield 3440.7 3426.0 14.7 




















Hot Hydrogen Materials and Component Tester 
Design Analysis - final selection 
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Hot Hydrogen Materials and Component Tester 
Conjugate Heat Transfer Analysis 
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Hot Hydrogen Materials and Component Tester 
Conjugate Heat Transfer Analysis - UNIC & SINDA Comparison 
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